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ABSTRACT 


In this research, we apply the time independent Schroedinger equation for a 
particle moving in one dimensional potential barrier of finite width and height 
We study the two cases which corresponds to the particle energies being 
respectively larger and smaller than the potential barrier. Then, we calculate 
transmission coefficient ( T) as a function of particle energy (£) for a potential 
barrier by changing the barrier height (Vo) and width (L] using Propagation 
Matrix Method. If we keep the barrier width constant and varying the height, 
we see that the passing limit is shifting towards the higher energies when 
barrier height is increased. If we keep the barrier height constant and change 
the barrier width, we see significance change in oscillations. 

KEYWORDS: Rectangular Potential Barrier , Propagation Matrix Method 

INTRODUCTION 

A differential equation for the wave function actually for the variation of 
wave function with space and time is called Schroedinger equation which 
describes the behavior of particles like electron, proton, neutron etc. We apply 
the Schroedinger equation for a particle moving in one dimensional potential 
barrier of finite width and height Vo. We consider a beam of particles of mass m 
along the x-axis from left to right on a potential barrier. According to classical 
physics, if a beam of particles with energy E < Vo is incident on the potential 
barrier, it will be reflected. It cannot go through the potential barrier[1]. 
However, according to quantum mechanics there is finite probability that a 
particle with energy less than the height of potential barrier can penetrate it. 


Tunneling phenomena are common atthe microscopic scale; 
they occur within nuclei, within atoms and within solids. In 
nuclear physics, for instance, there are nuclei that decay into 
an a particle and daughter nucleus[4J. The barrier 
penetration effect has important applications in various 
branches of modern physics ranging from particle and 
nuclear physics to semiconductor devices. For instance, 
radioactive decays and charge transport in electronic devices 
are typical examples of the tunneling effect[5J. 


The Rectangular Potential Barrier 

We consider a one dimensional potential barrier of finite 
width and height. The potential energy V (x) given by Eq.(l) 
is called the potential barrier which has a height of V 0 and a 


width of 

V(x) = < 


L. 

0 

Vo 

0 


x< 0 
0 <x<L 
x> L 


( 1 ) 


We consider particle of mass m incident on the barrier from 
the left with energy E. As mentioned therein, according to 
classical mechanics, the particle would be reflected back if E 
< Vo and would always be transmitted if E > Vo. We will show 
that, quantum mechanically, both reflection and 


transmission occur with finite probability for all values of E 
except in some special cases. We consider the following two 
cases which correspond to the particle energies being 
respectively larger and smaller than the potential barrier. 


A. Case I (£>14?) 

We divide the whole space into three regions: Region I (x < 
0 ), Region II (0 <x< L) and Region III (x> L). In region I and 
III the particle is free[2J. According to classical physics, if a 
beam of particles with energy E > V o approaches the 
potential barrier from the left, all of the particles in the beam 
will go over the barrier to region III.. The time independent 
Schroedinger equation for each region is 


— + k 2 
ax 


c) = 0 ( x < 0) 


( 2 ) 


' d 2 

K dx 2 



\f/ 2 {x) = 0 


[0 <x < L] 


' d 2 
v dx 2 



7 


^3 00 = 0 


(x>L) 


(3] 

(4] 


The solutions of the Schroedinger equation in the three 
regions are 
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y/(x) 


y/ l (x) = Ae ik ' x +Be~ ik ' x 

x<0 

y/ 2 (x) = C e iklX + D e~ iklX 

0 <x<L 

\f/,{x) = Fe k ' x 

x>L 


where k x = ^hnE/rf and k 2 =j2m(E-V 0 )/r] 2 .The 

potential barrier and propagation directions of the incident, 
reflected and transmitted waves are shown in Fig.(l) for E > 
Vo . The wave function will display an oscillatory pattern in 
all three regions; its amplitude reduces every time the 
particle enters a new region. The constants B, C, D and Fean 
be obtained in terms of 71 from the boundary conditions. The 
wave functions and their first derivatives must have 
continuous values [3]. Solving for F, we obtain 

F=4k l k 2 Ae lk ' L \(k l +k 1 ) 2 e~ lklL -{k^ -k 2 ) 2 e lklL Y 



Fig.l the potential barrier and the energy E 

F = 4 k x k 2 A e~ lklL \4k l k 2 cos(k 2 L)-2i(k x + k 2 )sm(k 2 L)\ l 

( 6 ) 


The probability of transmission is given by the transmission 
coefficient T. 

-i-i 


T = 


\A 


i+I 

4 


^ k\ -kl ^ 

k k 

v a 1 /v 2 J 


sin 2 {k 2 L) 


(■ E>V 0 ) (7) 


B. Case II (E < VO) 

According to classical physics, every particle that arrives at 
the barrier (x = 0) will be reflected back; no particle can 
penetrate the barrier. However, the quantum mechanical 
predictions differ sharply from their classical counterparts, 
for the wave functions is not zero beyond the barrier. In 
region I (x < 0 ), and region III (x > L), the Schroedinger 
equation and its solution remain the same as in case I. In 
region II [0 <x < L) the Schroedinger equation is 


' d 2 
K dx 2 



y/ 2 {x) = 0 


( 8 ) 


The solutions of the Schroedinger equation in the three 
regions are 

' ¥] (x) = Ae ik ' x + Be- lk ' x x<0 
y/(x) = \y/ 2 {x) = Ce klX + De~ klX 0 <x<L (9) 
y/ 3 (x) = Fe ik ' x x>L 


where k 2 =2mE/r\ 2 and k 2 = 2m(V 0 - E)/r\ 2 . The 

potential barrier and propagation directions of the incident, 
reflected and transmitted waves are shown in Fig.(2) for E < 
Vo. The wave function has an exponential form in the 
forbidden region inside the barrier. But there is also an 
oscillator wave to the right of barrier. To find transmission 


coefficients ^ _ 



we need only to calculate Fin terms 


of A. The wave functions and their first derivatives must 
have continuous values. Solving equation for F/A, we obtain 


E = 2e~ ik ' L 

A 


j2 _ 7 2 

2cosh(F 2 F) + z —- L sinh ( k 2 L ) 


k x k 2 


( 10 ) 



Fig (2) The potential barrier and the energy E . 


i i2 ^ 

4 cosh 2 (k 2 L ) + 

( kl-kF 

2 

sinh 2 (k 2 L) 

\ A \ 


{ k Ai ) 



( 11 ) 


Since cosh 2 (k 2 L) = 1 + sinh 2 (£ 2 L) , we can reduce to 


"i + i 

fk 2 + kF 

4 

k k 

V ^2 J 


sinh 2 (k 2 L) 


( 12 ) 


We note that T is finite. This means that the probability for 
the transmission of the particles into the region 0 < x < L 
is not zero as expected from classical physics. This is a purely 
quantum mechanical effect which is due to the wave aspect 
of microscopic objects; it is known as the tunneling effect: 
quantum mechanical objects can tunnel through classically 
impenetrable barriers. 

The Propagation Matrix Method 
The Propagation Matrix 

A method is needed a method for finding solutions to 
complicated potential structures for which analytic 
expressions are unmanageable. The transmission coefficient 
is calculated at the first potential step for a particle of energy 
E incident from the left. We then imagines the transmitted 
particle propagating to the next potential step, where it 
again has a probability of being transmitted or reflected. 
Associated with every potential step and free propagation 
region to the next potential step is a 2x2 matrix which 
carries wave function amplitude. The total one dimensional 
propagation probability for a potential consisting of a 
number of potential steps may be calculated by multiplying 
together each 2x2 matrix associated with transmission and 
reflection at each potential step. Therefore, the wave 
function coefficients for a particle traversing a one 
dimensional potential consisting of a number of such regions 


@ IJTSRD | Unique Paper ID - IJTSRD26813 | Volume - 3 | Issue - 5 | July - August 2019 


Page 1892 
















































International Journal of Trend in Scientific Research and Development (IJTSRD] @ www.ijtsrd.com elSSN: 2456-6470 

may be calculated by multiplying together the appropriate 2x2 matrices. We can solve for a particle moving in an arbitrary 
potential by dividing the potential into a number of potential energy steps. The following four basic parts are needed. We may 
use the propagation matrix method to calculate the probability of the electron emerging on the right-hand side of the barrier. 
The method is best approached by dividing it into small, easy-to-understand, logical parts. 


The Step Propagation Matrix 

We calculate the propagation matrix p s tep for transmission and reflection of the wave function representing a particle of energy 
E incident on a single potential step. The potential step we consider is at position Xj+i in Fig.(3]. 



Position, x 

Fig.3 Approximation of a smoothly varying one dimensional potential V(x] with a series of potential steps 

In this approach, the potential between position Xj and Xj+i in region j is approximated by a value Vj. Associated with the 
potential step at xj and free propagation distance Lj = Xj+i - Xj is a 2x2 matrix which carries all of the amplitude and phase 
information on the particle. Fig.(4] shows detail of the potential step at position index j +1. The coefficients A and C correspond 
to waves traveling left to right in region j and j +1, respectively. 


A 



n 


Energy j E 


Regionj 


Regionj + 1 


V i + i 


Position, x x_ x j + 1 

Fig.4 A one dimensional potential step. In region j the potential energy is Vj and in region j+1 the potential energy is Vj+i. 


The transition between 

„ _ ( 2 m (E - Vj)) 1 ' 2 

R 

j+1, are 


region 1 and region 2 occurs at position x = Xj +i. The electron (or particle] has wave vector 
in region j, and the wave functions, which are solutions to the Schrodinger equation in regions j and 


Vi= A ,e 


ik, 


+ 


b j 


-i ki 


(13) 


¥ 


7 + 1 


= c Jtl 


i k 


7 + 1 - 


+i v. 


-i k 


j + 1 - 


(14) 


Following the convention we have adopted in this paper, A and C are coefficients for the wave function traveling left-to-right in 
regions j and j+1, respectively, and B and D are the corresponding right-to-left traveling-wave coefficients. The two wave 
functions given by Eq.(13) and (14] are related to each other by the constraint that i |j and d^/dx must be continuous. This 
means that at the potential step that occurs at the boundary between regions j and j+1 we require 
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and 


A y+1 

dy, 


= ¥ j+ 1 


d x 


x=x j + i 

d¥, + i 


d x 


A j +1 


CIS) 

(16) 


•*y + i 


Substituting Eq.(13] and (14] into Eq.(15] and (16] gives two equations 

A ik ; x t~j -ik;X ^ ik ; ,iX ~ik:. iX 

V ' +B ; e ' =C Jtl e '*' + £>, +1 e 


a ik;X ^ , ,v, ^ 

Aj e J -Bj e J - 


ikjX _ AV y + l 


_C g ^ 

k; i + l 


k 

1 k j+ ]X 7+1 


-/^• +1 x 


_n ^ ^ 

k, ;+i 


(17) 

(18) 


By organizing into rows and columns the terms that contain left to right traveling waves of the form e ikx and right-to-left 
traveling waves of the form e~ ikx , we may write Eq.(17] and (18] for a potential step at position Xj+i= 0 as a matrix equation: 
1 1 

k J+ 1 k J+ 1 ~ ;+1 (19) 


1 1 
1 -1 


A. 


B , 


k. 




-1 



We would much prefer a simple equation of the type 


~ A J~ 


~c J+ T 

B 

P j step 

D fx . 

L 7 J 


L 7 + 1 J 


( 20 ) 


where Pj S te P is the 2x2 matrix describing wave propagation at a potential step. To obtain this expression, we need to eliminate 
the 2x2 matrix on the left-hand side of Eq.(19]. We simply use from basic linear algebra that the inverse of a 2x2 matrix 


A = 


a lx 

a \ 2 

is 

A _1 


l 

a 2 2 

a 1 2 

a 21 

a 22 

si 


| A 

— a 2 j 

a i i 


|^41 — CL xl 2 ~ a \ 2 l ■ Hence, the inverse of 

1 







1 

"1 1 " 


5 ,. 

2 

1 - 1 _ 



1 

k j +1 


j + i 
kj 


C J + i 


D 


j + i 


step 


where the determinant of A is given by 

, so that we may write 

( 21 ] 


“1 

1 “ 

is 

1 

- r 

i 

"1 

1 " 

1 

-1_ 

2 

1 

i 

2 

1 

- 1 


j + i 

j +1 


where the step matrix is 

k . 


1 


j step 


1 + 


k ) 


1 - 


k j 




k j 


k , 


( 22 ) 


This is our result for the step potential that will be used later. We continue the development of the matrix method by 
considering the propagation between steps. 

The Propagation between Steps 

We calculate the propagation matrix pf ree for propagation of the wave function between steps. The free propagation we 
consider is between positions Xj and Xj+i in Fig.(3]. The distance of this free propagation is Lj. Propagation between potential 

steps separated by distance Lj carries phase information only so that A .e J ’ =C + 1 andi? e lk] 1 = Z ) +1 .This maybe 


expressed in matrix form as 


i kj L: 

e 1 1 

i- 

o 



1 

+ 

_ 1 

i 

o 

-i k-L ■ 

e 1 1 

_ b j_ 


l D J^\ 


Or, alternatively, 
where p jfree = 


[A.l 


rc, + ii 

j 

B J_ 

B j free 

j +1 

_^7 + l_ 


(23) 

(24) 

(25) 
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The Propagation Matrix Pj for they-th Region 

We calculate the propagation matrix for the j-th region in Fig.(3). This is achieved if we multiply p s tep and pf ree to obtain the 
propagation matrix Pj for the j-th region of the discretized potential. To find the combined effect of pf ree and pj S tep we simply 
multiply the two matrices together. Hence, propagation across the complete j-th element consisting of a free propagation region 
and a step is 


Pj = Pjfree P. 


jstep 


Pl 1 

Pl 2 

Pl 1 

Pll 


(26) 


When we multiply out the matrices Pjfree Pjstep given by Eq.(25) and (22), respectively, it gives us the propagation matrix for the 
j-th region: 


i + ^ 

v k J 

k ^ 

i _ zm 

k j j 


f k ' 

i _ */±i 

k j j 

k. ^ 

k j j 


P 11 

Pn 

Pn 

* 

Pi 1 


(27) 


Propagation through an Arbitrary Series of Step Potentials 

We calculate the total propagation matrix P for the complete potential by multiplying together the propagation matrices for 
each region of the discretized potential. For the general case of N potential steps, we write down the propagation matrix for 
each region and multiply out to obtain the total propagation matrix, 

P = Pl p 2 ...pj ...p N = fj \ Pj (28) 

7 = 1 


The total propagation matrix P satisfies continuity in *¥ and dW /dx between adjacent regions. Since the particle is introduced 
from the left, we know that A = 1, and if there is no reflection at the far right then D = 0. We may then rewrite 


~ A " 


( j = N A 

“ c“ 


"C “ 

B 

— 

TlPj 

W = 1 J 

D 

= P 

D 


as 

" 1 " 


" Pn 

P 12 ’ 


" c" 


B 


P 21 

P 22 


0 


(30) 


In this case, because 1 = pn C, the transmission probability I C | 2 is simply 


C = 


1 


(31) 


Eq.(20) is a particularly simple result. We will make use of this when we calculate the transmission probability of a particle 
through an essentially arbitrary one-dimensional potential. 


Transmission Probability for a Rectangular Potential Barrier 

Fig.(5) is a sketch of the rectangular potential barrier we will consider. The thickness of the barrier is L. A particle of mass m 
incident from the left of energy E has wave vector ki. In the barrier region, the wave vector is k 2 . The wave vector ki and k 2 are 

related through k\ — k\ + 2 m V 0 / T| 2 . 



£P 

Q 

O 


O 

o 



X = Q X = L Position, X 

Fig.5 The potential of a one dimensional rectangular barrier of energy Vo 


@ IJTSRD | Unique Paper ID - IJTSRD26813 | Volume - 3 | Issue - 5 | July - August 2019 


Page 1895 













































International Journal of Trend in Scientific Research and Development (IJTSRD] @ www.ijtsrd.com elSSN: 2456-6470 

A wavy particle incident on the barrier from the left with amplitude A sees a potential step-up in energy of V 0 at X = 0 , a 
barrier propagation region of length L , and a potential step-down at X = L . A particle of energy E , mass m , and charge e 

has wave number k x outside the barrier and k 2 in the barrier region 0 < X < L .We consider a particle impinging on a step 

change in potential between two regions in which the wave vector changes from ki to k 2 due to the potential step up shown in 
Fig.(3). The corresponding wave function changes from i|/i to i|/ 2 . Solutions of the Schroedinger equation for a step change in 
potential are 


^i = 


w 2 = 


A 

Eh 


i ki x 

e 1 + 


B 

Eh 

D 

Eh Eh 


C i 


+ ■ 


-i k x x 


-i k 2 x 


Applying the condition that the wave function is continuous at the potential step W x \ step — *P 

dW dE 

of the wave function is continuous-- =-- gives 


\step 


(32) 

(33) 

and that the derivative 


step 


step 


A B 
- + ■ 


C D 
■ + ■ 


"E\ ■\fk y[k 2 " 

A B k , C k , 


D 


Rewritten in matrix form, these equations become 


(34) 

(35) 


1 

'l 

1 ' 

~A 

1 

" i i 

k „ k o 

"c" 


_1 

- 1 _ 

_B_ 

Ei 

-1 

s |^“ 

l 

l^r 

_i 

D 


(36) 


To eliminate the 2x2 matrix on the left-hand side of this equation, we must find and multiply by its inverse matrix. The 
determinant ofthe left-hand matrix is (— 1 — \)/h = — 2/ k x , so the inverse of the left-hand matrix is !kl ^ ' .Hence, 

2 VU 1 ~ l . 

we may rewrite Eq. (36) as 


' A ' 

k x 1 

'l l ' 

' l l 

k, k, 

1 

"c" 

1 

'i i' 

k ! k x 

"C" 

_B_ 

(N 

1 

! -l_ 

1 - 

| n 

1 

1_ 

E 

D 

2 yjk x k 2 

1 -i_ 

k 2 - k 2 

D 


Multiplying out the two square matrices gives the 2x2 matrix describing propagation at the step-up in potential 


B 


1 


k x + k 2 


k x - k 2 


2jkA~lk l -k 2 k x + k 2 


C 

D 


(37) 


Since the rectangular potential barrier consists of a step up and a step down, we can make use of this symmetry and 
immediately calculate the 2x2 matrix for the step down by simply interchanging ki and k 2 . The total propagation matrix for 
the rectangular potential barrier of thickness L consists of the step-up 2x2 matrix multiplied by the propagation matrix from 
the barrier thickness L multiplied by the step-down matrix. Hence, our propagation matrix become 

p = 


1 

K + k 2 k \~ k 2 

-1 

tM 

0 

1 

k 2 +k x k 2 - k 

2 El k 2 

k x - k 2 k x + k 2 

0 

1_ 

2-yJk\ ^2 

k 2 ~k k 2 + k x 


p = 


1 

\k+k 2 )e~ ihL (k-k 2 )e n > L ~ 

k 2 +k x k 2 ~k 

4 k { k 2 

(k - k 2 )e~ l k - L (k +k 2 )e lklL 

k 2 ~k k 2 + k 


(38) 


To find the matrix elements of P, we just multiply out the matrices in Eq.(38). For example p n becomes 


k 2 -k 2 

Pn = - 2 1 ( e ik ' L -e~ ihL ) 

12 4 k x k 2 


(39) 
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The next step we want to take is to calculate the transmission probability for a particle incident on the barrier. We already 


know that the transmission of a particle incident from the left is given by 
obtaining p n from Eq.(38) 


1 


, so that we will be interested in 


Pi i = 

Pi i = 

Pn = 


(* 2 +*i)(*i +kpe^+(k 1 -k^jk, M 

4k x k 2 

(k 2 + k x + 2k x k 2 )e~ lklL —(k x +k 2 —2k x k 2 )e lklL 

4k x k 2 

(kl + k x )(e~ iklL — e iklL ) | 2 k x k 2 (e~ ik2L + e iklL ) 
4/c, k 2 4 k x k 2 


2 2 k x k 2 

2 2 k x k 2 


^{e~ ik > L +e ik > L ) 

-( e - ik > L +e ihL ) 

2 


(40) 

(41) 

(42) 

(43) 

(44) 


The Result and Discussion 

We use the propagation matrix approach to calculate the 
transmission coefficient T of a rectangular potential barrier 
with width L = lnm and height Vo = 0.3 eV as shown in 
Fig.(5). Matrix technique in transmission calculations means 
that we express every region of barrier as a propagation 
matrix P. The transmission coefficient can be expressed 


T 


1 . Using this principle, transmission coefficient T 

i 2 

P 11 | 


dependence on electron energy for a rectangular potential 
barrier was numerically calculated in Fig.(6). We vary 
barrier width while keeping barrier height constant. 
Transmission probabilities T ( E) for barrier height Vo = 0.3 
eV and width L = 0.2 nm to 3.2 nm by increasing twice of 
initial value are shown in Fig.(7) to Fig.(ll) respectively. 
Fig.(12) shows the transmission probability T (. E ) for barrier 
height Vo = 0.3 eV and different barrier width. 


Transmission probability fcr a potential barrier of width 1 nm and height 0.3 eV 



Fig.(6) The transmission coefficient T for a barrier of 
width L = 1 nm and height Vo = 0.3 eV. 


Transmission probability for a potential barrier of height 0.3 eV and Length 0.2 nm 



Fig.(7) Numerical calculation of transmission coefficient 
for a barrier of height Vo = 0.3 eV and width L = 0.2 nm. 


Transmission probability for a potential barrier of height 0.3 eV and Length 0.4 nm 



Fig.(8) Numerical calculation of transmission coefficient 
for a barrier of height Vo = 0.3 eV and width L = 0.4 nm. 
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Transmission probability for a potential barrier of height 0.3 eV and Length 0.8 nm Transmssion probability for a potential barrier of height 0.3 eV and different width 



Fig.(9) Numerical calculation of transmission coefficient 
for a barrier of height Vo = 0.3 eV and width L = 0.8 nm. 


Transmission probability for a potential barrier of height 0.3 eV and Length 1.6 nm 



Fig.(10) Numerical calculation of transmission coefficient 
for a barrier of height Vo = 0.3 eV and width L = 1.6 nm. 


Transmission probability for a potential barrier of height 0.3 eV and Length 3.2 nm 



Energy [eV] 


Fig.fll) Numerical calculation of transmission coefficient 
for a barrier of height Vo = 0.3 eV and width L = 3.2 nm. 



Fig.(12) Numerical calculation of transmission coefficient 
for a barrier of height Vo = 0.3 eV and different width. 


When a beam of particles of fixed energy is incident on a 
potential barrier, a certain fraction of the incident particles is 
transmitted while the remaining fraction is reflected. This is 
in contrast to classical mechanics. Classically, there must be 
total transmission if the energy of the incident particle is 
more than the height of the barrier, and total reflection if the 
energy of the incident particle is less than the height of the 
barrier. Fig.(12) shows the transmission probability T[E) for 
barrier height Vo = 0.3 eV and different barrier width. If we 
keep the barrier height constant and change the barrier 
width, we see significance change in oscillations. When 
barrier width is narrow, even electron with lower energy can 
pass through the barrier by quantum tunneling. The wider 
the barrier will be the less observable the tunneling is, which 
is expectation from classical approach as well. 

Conclusion 

We have found that the transmission coefficient [T] depends 
upon potential barrier height [Vo] and barrier width ( L). The 
transmission coefficient (7) is a measure of the probability 
that the particle will be transmitted through the barrier. 
Thus, we conclude that there is finite probability of particle 
penetrating the barrier and appearing on the other sides. 
The ability of particle to penetrate the barrier when E < Vo is 
a quantum mechanical result and is known as tunnel effect. 
In nuclear physics, there are nuclei that decay into an a 
particle and daughter nucleus. This a decay process can be 
explained by tunnel effect. 
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